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114 cross sections for nuclide production in a 1.0 GeV proton-irradiated thin ^"^Pb target have been 
measured by the direct 7-spectronietry method using a high-resolution Ge detector. The 7 spectra 
were processed by the GENIE-2000 code. The ITEP-developed SIGMA code was used together with 
the PCNUDAT nuclear decay database to identify the 7 lines and to determine the cross sections. 
The ■^^Al(p,x)'^^Na reaction was used to monitor the proton flux. Results of a feasibility study of 
the auxiliary ^'^AI(p,x)'^'*Na and ^'^Al(p,x)^Be monitor reactions in the 0.07-2.6 GeV proton-energy 
range are presented as well. Most of the experimental data have been analyzed by the LAHET 
(with ISABEL and Bertini options), CEM95, CEM2k, INUCL, CASCADE, CASCADE/INPE, and 
YIELDX codes that simulate hadron-nucleus interactions. 
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I. INTRODUCTION 

Major advances in Pb-Bi reactor technology, as mani- 
fested in the Russian-made Alpha-class submarines, has 
made Pb-Bi technology attractive for accelerator-driven 
systems (ADS) to be used at facilities designed to trans- 
mute nuclear waste These Pb-Bi reactors would 
operate in an environment of high-energy radiation, 
thus motivating studies of the nucleonics characteristics 
of lead and bismuth, including the yields of residual 
product nuclei under irradiation by protons of energies 
ranging from a few MeV to 2-3 GeV. Results of such 
studies are extremely important when designing even 
demonstration versions of ADS facilities. Yields of 
residual product nuclei define important characteristics 
of these facilities such as radioactivity (both current and 
residual), deterioration of corrosion resistance, yields of 
gaseous products, neutron "poisoning", etc. That is why 
several scientific groups have been studying the yields of 
residual product nuclei from Pb and Bi (both natural and 
isotope-enriched) bombarded with protons in the desired 
energy range |Q-|lO|. 

Most of these studies have been made using 7 spec- 
trometry on proton-irradiated natural or monoisotopic 
thin Pb and Bi samples (for example, p-l-^glPb). This 
technique makes it possible to determine the independent 



and cumulative yields of more than a hundred reaction 
products An alternative technique of magnet- 

ically separating the reaction products after a heavy- 
ion beam interacts with a liquid H target, using reverse 
kinematics, as done recently at GSI to study ^gfU, ^ygAu 
and ^g|Pb-|-}H, makes it possible to determine about a 
thousand yields 

Comparing the results of the two techniques helps in 
understanding the systematic errors inherent to both. 
This work makes use of the 7-spectrometry technique. 
We discuss below our results as well as the systematic 
errors of 7 spectrometry. 

Necessarily, computational methods will play an im- 
portant role in the designs for ADS facilities. Therefore, 
testing of the most extensively used simulation codes is 
necessary. We perform a qualitative and quantitative 
comparison of our data with simulated results from seven 
codes widely used in applications. 



II. BASIC DEFINITIONS AND FORMULAS 

The formalism for finding the reaction product yields 
(cross sections) in high-energy proton-irradiated thin 
targets is described in sufficient detail in Any of 
the measured reaction products generated in nuclcon 
interactions with matter are assumed to originate both 
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in the reaction proper and in the decays of its chain 
precursors. Thus, the set of terms invented earher in 
studying the mass and charge distributions of fission 
products (see, e. g., ijlll) can conveniently be used when 
processing and interpreting the experimental results of 
the present work. In this terminology, the independent 
and cumulative yields of products underlie the formalism. 

The independent yield of a reaction product nuclide is 
a cross section for the nuclide to be produced directly in 
the reaction, whereas the cumulative yield of the nuclide 
is a cross section for the nuclide to be produced in all the 
appropriate processes, i.e., both directly in the reaction 
and over the time during the decays of all of its chain 
precursors. [] 

The variations in the concentration of any two nuclides 

of a chain produced in an irradiated target (A^i 

N2 may be modeled as a set of differential equations 
that describe the production and decays of the nuclides. 
By introducing the time functions, Fi, of the type Fi = 

(1 — e~^''^) ^ _x T ' (* = I7 2, or monitor product 

(Na or another)); r is the duration of a single proton 
pulse, T is the pulse repetition period, K is the number of 
irradiation pulses which characterize the nuclide decays 
within the irradiation time, and by expressing (similar 
to the relative measurements) the proton fluence via the 
monitor reaction cross section agt, we can present the 
cumulative and independent yields as 



cum 



Aq NaI F^a 

rjiEiFiNNa Nt XNa 



(1) 



or, if measurements of the second nuclide alone are used, 
as 



And 



Ai NaI A2 ~ Ai FNa 

VimSlFlNNa Nt A2 \Na 
M _|_ ^Ai\ 1 Nai FNa 

F2 Fi X2 J ri2^2NNa Nt XNa 
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Nai FNa 



Fi F2 J ri2e2NNa Nt XNa 



(2) 
(3) 



(4) 



Here, tr™™ is the cumulative cross section of the first 
nuclide iVi; a™'^ and ctj"™ are the independent and 
cumulative cross sections of the second nuclide N2] Nai 



*Some applications also make use of the mass yield. In its 
simplified form (disregarding alpha transitions and delayed- 
neutron emission), the mass yield is the sum of all indepen- 
dent yields of the isobars of a given mass or equivalently is 
the sum of cumulative yields of the stable isotopes of a given 
mass. 

^In case of absence of short-lived precursors, the yield of the 
nuclide A'^i should be considered as independent. 



and Nt are the numbers of nuclei in the monitor and in 
the experimental target, respectively; NNa is the number 
of Na nuclei produced in the monitor; r/i and 772 are 7- 
line yields; ei and £2 are the spectrometer efficiencies 
at energies Ej-^^ and Ej^ ; vi is the branching ratio of 
the first nuclide; Ai, A2, and XNa are, respectively, the 
decay constants of the first and second nuclides and of 
the monitor product (-^^Na and/or ^"'Na). 

The factors Aq, Ai, and A2 are calculated through 
fitting the measured counting rates in the total absorp- 
tion peaks, which correspond to energies Ej-^ (the first 
nuclide) and E-y^ (the second nuclide), by the exponential 
functions 



git) = ^0 



1 - e-^i**-- 
Xitt rue 



-Alt 



(5) 



and 



fit) = A: 



A, 



1 - e-^i*'™ 

Xltfrue 



-Alt 



X2ttr 



(6) 



Here, ttme is the real time over which the spectrum is 
measured, t is the time between the end of the irradiation 
and the beginning of the spectrum measurement. The 
forms of functions (|^) and (^ are the similar to what is 
presented in but we present them here for clarity: 

Ao = 7VT$^Tr™ryieiFi , (7) 
Ai = NT^a1''"'m£2Fii^iX2/iX2 - Ai) , (8) 
A2 = NT<^mF2K''' - <TT"'i^iXi/iX2 - Ai)] . (9) 

It should be noted that formulas (|l|)-(0) are deduced 
on the assumption that the corresponding 7 counting 
rates of each nuclide produced under irradiation are 
determined to within the desired accuracy throughout 
the time interval from the irradiation end to the ultimate 
intensity detection threshold. Curve 1 of Fig. |l| exempli- 
fies such a favorable situation. The situation corresponds 
to the i^^Hg (Ti/2 = 4.85 hours) ^ ^'^'^ku (T1/2 = 4.94 
hours) decay chain, with the ^^^Au, 316.5 keV 7 line 
being measured. Despite the very similar periods, the 
absence of any addends to the measured 7 line (i.e., null 
contribution from 7 lines of any other nuclides whose 
energy is the same, within the spectrometer resolution, 
as of the measured 7 line) provides a sufficiently accurate 
determination of the cumulative ^^^Hg yield, as well as 
the independent and cumulative ^^^Au yields. 

The real situation often gets more complicated, how- 
ever, because the decay curve of the second nuclide can- 
not be measured correctly in every case within the desired 
(in the sense as mentioned above) time interval. Its 7-line 
intensity is often difficult to measure in the beginning of 
a measurement run, right after the irradiation stops. In 
the case of short-lived nuclides, this is due to the fact 
that within the period from the end of irradiation to 
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(1) E =316.5keV ^^^Hg(4.85h)^^^^Au(4.94h) 

(2) E^=2214.6keV ^^^Pt(10.2d)^^^^lr(41.5h) 

(3) E^=123.7keV ^^^Ta(3.14h)^^^^Hf(23.6h) 

(4) E =160.4keV^^^Ta(3.14h) 

(5) E'=172.2keV ^^^Ta(3.14h) + ^^>t(2.9d) 




2 3 4 5 6 
Cooling time [days] 

FIG. 1. Typical examples of decay curves. Curve (1) 
is for the chain ^^^Hg -> ^^^Au. Curve (2) is for ^^^Pt ^ 
i**Ir. Curve (3) is for ^''^Ta ^ "^'Hf. Curve (4) is for the 
independent ^^"^Ta decay. Curve (5) is for the independent 
(^^^Ta + ^^^Pt) decay. Curve (1) is drawn with a time scale 
factor of 4 and a counting rate factor of 0.1. Curve (2) is 
drawn with a time scale factor of 0.1 and a counting rate 
factor of 100. 



the beginning of a measurement run (the coohng time), 
the nuchde A^i (precursor) can partly or fully decay (if 
Ai > A2), or else a full or partial equilibrium sets in (if 
Ai < A2) between the nucHdes A^i and N2. 

The contribution from nuclide Ni will then never 
be reflected in the experimental decay curve of nuclide 
N2- The situation gets even more complicated because 
a great number of the products can include nuclides 
whose half-lives are close to the half-life of a shorter- 
lived nuclide, either iVi or N2- In that CcLSG, clS noted 
above, if the 7-line energy of such nuclides is the same 
(within the spectrometer resolution) as the measured 7- 
line energy of other nuclide, then neither of the factors 
Ai and A2 can be determined, or else the two can only 
be determined with great uncertainty and thus become 
useless for calculating the yields. 

Fig. |l] is also a good illustration of a possible unfa- 



vorable situation for analyzing the nuclides of similar 
half-lives (Curve 1). If the ^^^Au decay curve begins 
being measured more than two days after the irradiation, 
the ^^^Hg contribution becomes uncertain, resulting in 
erroneous calculations of the ^^^Au yield. Such a 
situation is most probably responsible for the difference 
in the measured ^^^Au yields of more than a factor 3 
between the data of the present work (see Table II below) 
and the results of g] (46.9±6.6 mb and 160±50 mb, 
respectively) . 

Analyzing possible structures of radioactive chains 
permits the following two very common situations to be 
singled out. 

First, assume that Ai < A2. This situation is exempli- 
fied in Fig. |l| (Curve 2), which shows the decay curve of 
chain nuclides i88pt(Ti/2=10.2 days)^i8Slr(Ti/2=41.5 
hours) recorded by measuring the 2214.6 keV 7 line of 
the is^Ir daughter nuclide. In this case, the decay curve 
of nuclide i^^Ir can be used to obtain fairly accurate 
values of the factors Ai and A2 and, hence, of crj" 



And 



and (To 



cum 

88 



-ind 

' 188/^1 



and af|'s7J. Should 



this favorable situation get complicated (for example, 
the measurements began in 2 days after irradiation), 
then, even without observing the knee, the conclusion 
concerning the ^®*Pt production is quite obvious because 
the 2214.6 keV 7 line of ^^^Ir is measured with the ^^^Pt 
period. In this case, formula (0) may be used to calculate 
the crf"s™j value, whereas the Ir data prove to be lost. 

It should be noted that there even exists a nuclear data 
library that presents the ^*^Pt gamma yields corrected 
for the (A2 — Ai)/A2 value. The corrected yields are 
marked with an index (D) to notify the validity of using 
the daughter nuclide 7 lines when calculating the number 
of mother nuclei ||l^ . 

The inverse situation (Ai > A2) permits the factor 
A2 alone to be determined reliably. This is exempli- 
fied by Curve 3 in Fig. 0, which is the decay curve 
of the chain nuclides ^^^Ta (Ti/2=3.14 hours)^"3jjf 
(Ti/2=23.6 hours) recorded by measuring the 123.7 keV 
7 line of ^^•^Hf. In this case, the factors A\ and 
A2 can be determined to within the desired accuracy, 
and formulas (||)-(0) can be used to calculate crj"™, 
and cr^"" (CTf7"3Ta' ^i?^'^/' and afr'^j'^.). In cases 
when unfavorable situations arise due to, for example, 
measurements starting a day after the irradiation, the 
factor Ai alone can be determined. 

If, however, the eigen 7 line is used (this is shown 
in Fig. |l| by Curve 4, which is the ^^^Ta decay curve 
measured by the 160.4 keV -'^^^Ta 7 line, or by Curve 5, 
which is the same, but inferred from the 172.3 keV ^^'^Ta 
7 line with an addend contributed by the ^^^Pt decay), 
then the factor Aq can be calculated and the missing 
factor Ai is 



Ai = Ao vi 



A2 



A2 — Ai 

can be used to find crj 



(10) 



whereupon formulas (Q 

a|"^ and ct^"" (<^™3Ta' ^"3"^/. and af^"^^.). If, however. 
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the factor Ai cannot be found, we may use the factor A2 
together with expression (14) from ||6| to determine the 
quantity cr™™ , which we call the supra cumulative yield: 



ind 



+ 



Ai 



cum 
Via I 



Ai — A2 

A2 NaI FNa ^ 



(11) 



From the resultant formula (|T^) it is seen that the 
measured value of the supra cumulative yield cr™"* may 
sometimes prove to be very different from its true value 
fjcum^ In the case of ^^^Re, for example, the data of 
the present work give a ^ctj"™ value of ^ 55%. This 
fact should be borne in mind when comparing between 
experimental and simulated data. 



It should be noted that the difference between tr™™ 



and 



is not specified in many of the relevant 



publications despite the fact that ctj"™ is always greater 
than (Tj"™- The explanation is that the formal equality 
g-cum _ ^cum j^Q^g jj^ ^g^gg |;]^g generation times of the 

first and second nuclides are different. In this case, the 
time correction can be determined: 



A, 



A. 



(12) 



From formula (|12|), it is seen that At depends on 
the yields erf"™ and cr^'^j thereby preventing the time 
correction Ai from being allowed for in a general case 
when determining 0-2"'". In common cases of cr™'' ^ 
i/icrj"'", however, the time correction At can be shown 
to depend on the decay constants only: 



At ^s: At' = -!- 
A2 



1 



A, 



Ai — A2 



(13) 



Therefore, the cumulative yield of the second nuclide 
can be determined after a post-irradiation period suffi- 
cient for A^i to decay into N2 (normally, this equals from 
6 to 10 half-lives of the first nuclide) by measuring the 
decay curve of the first nuclide making allowance for the 
time correction At: 



A. 



ri2S2F2NNa Nt XNa 
A2 NaI FNa 

r]2S2F2NNa Nt XNa 



NAlFNa -X^At 



CTst 1 



A2 

Ai 



(14) 



In the case where the condition cr™ <^ (t™"Vi is not 
satisfied, so the cr™™ value cannot be calculated in any 
way accurately, we may estimate the difference 



AaS 



Ao 



Ai-A: 



cum 
-1^1 (J I 



(15) 



Proceeding from the condition cr™' 

we may estimate the upper limit of Acrf 



> erf"" > 



crf""j/i 



III. EXPERIMENTAL TECHNIQUES 

The used experimental techniques are described in [|| . 
A 10.5 mm diameter, 139.4 mg/cm^ monoisotopic ^"^^Pb 
metal foil sample (97.2% ^ospb, 1.93% 207pb, 0.87% 
206pb, <0.01% 204pb, < 0.00105% of chemical impuri- 
ties) was irradiated by protons. A 139.6 mg/cm^ AI foil 
of the same diameter was used as monitor. Chemical 
impurities in the monitor do not exceed 0.001%. 

The sample was irradiated by the 995 MeV external 
proton beam of the ITEP U-10 synchrotron described 
briefly in |^ . The irradiation time was an hour. The total 
proton fluence was 5 • lO^'^p/cm^. The 7-spectrometer 
measurements began 10 minutes after the irradiation 
stopped. The spectra have been measured for half a year. 

It should be noted that the experimental program 
of measuring the reaction product yields from diverse 
materials was supported by other studies aimed at 
reducing the systematic errors in the experimental results 
p3[ . These studies include: 

• experiments to specify the neutron component in 
the extracted proton beams, 

• experiments to specify the monitor-reaction cross 
section, 

• experiments to specify the variations of the 7- 
spectrometer detection efficiency, depending on the 
height of the irradiated sample position above the 
detector, 

• studies to optimize the 7-spectrum processing 
codes. 

The results of these studies were used in determining 
the cross sections for nuclide production in the measured 
^*'^Pb(p,x) reactions. 



Neutron component in the extracted proton 
beams 



< 



Ai — A2 



or, in the relative form, 



Act™™ 

"2 



• 100% < 



Ai - A, 



• 100% 



(16) 



(17) 



The proton beams extracted from accelerators in- 
clude not only the primary protons, but also secondary 
particles (neutrons, protons, tt mesons, and gammas) 
produced in the primary proton interactions with the 
structural materials of transport channels and shielding. 

Identical reaction products can be produced in interac- 
tions of various secondaries with an experimental sample. 
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Since particular nuclear reactions that generate a given 
nuclide cannot be identified in the measurements, the 
extracted proton beams have to be tested and specified 
thoroughly. 

Track detectors were first used in the experiments 
to discriminate the neutron component in the proton 
beams Later, direct 7 spectrometry was used for the 
purpose. Solid-state nuclear track detectors (SSNTD) of 
an improved geometry with a coUimating grid and a glass 
fragment track detection material were used to record the 
fission fragments from a fissile layer, thereby improving 
the absolute detector efficiency. 

An SSNTD with a 61.5 /xg/cm^ ^°^Bi layer was used to 
measure the proton flux density. ^°^Bi was selected be- 
cause the cross section for its fission induced by secondary 
neutrons is small compared with that for primary protons 
{(7209 Bi{n.f) << o'209s,;(p,/)). The neutron flux density 
was measured using an SSNTD with a 880 /ig/cm^ ^'^^Np 
layer. 

The following experimental design was adopted. The 
extracted proton beam irradiates a ^°^Bi-containing 
"sandwich" (Bi layer -I- collimator -I- glass), while similar 
sandwiches with ^'^^Np layers are placed along a line 
normal to the beam axis at distances of 20-435 mm from 
the axis. 

In the experiments, the neutron-to-proton flux density 
ratio, <i>„/$p, was measured using the expression 



T2 



N 



•'Bi 



(18) 



where Ti and T2 are numbers of measured tracks of 
^^^Np and ^°^Bi fission products, respectively; N 
and A'^ are numbers of the ^'^^Np and ^°^Bi nuclei, 
respectively; ^1 and ^2 are, respectively, corrections to the 
^^^Np and ^"^Bi layers, which allow for the anisotropy 
of fission-fragment ejection and for the variations of 
the solid angle of fission-fragment ejection through the 
collimator grid; is the cross section for proton- 

induced ^"^Bi fission; ^ is the weighted mean ^'^'^Np 
neutron-induced fission cross section calculated as 



_ ^ / a^{E)^{E)dE 



(19) 



where 

X stands for 237Np(n,f), 27Al(n,p)27Mg, 27Al(n,x)24Na, 
and 27Al(n,x)22Na. 

The fission cross section cr^ ^ (E) needed to calcu- 

late (Tjj J was retrieved from the WIND data library 
The cross section for the proton-induced ^^^Bi 
fission, <7pj^\ was taken from | p^ . 

The experiments were made with 200, 800, and 2600 
MeV proton beams. Fig. |^ shows the resultant $n,/$p 
ratios as functions of the perpendicular distance to the 
proton beam. The $„/$p ratio right in the proton beam 



was estimated by extrapolating the peripheral results to 
the center and proved to be of about (0.2-2)%. 

The feasibility of distinguishing the (n,p) reactions 
from (p,x) reactions has permitted an alternative pat- 
tern of direct 7 spectrometry. Namely, Al sam- 
ples were irradiated, and 2^Al(n,p)2'''Mg (a ~ 2.5MeV 
threshold), p'^Al(n,a)24Na (a ~5.5 MeV threshold) -t- 



27Al(p,: 



\24 



Na 



(a -25 MeV threshold)], ^^A\{^,x 



\22 



Na, 



and 2^Al(p,x)^Be reaction rates were measured, in the 
beam center and at periphery. The 2^Al(n,p)2^Mg reac- 
tion characteristics have made it possible to detect ^^Mg 
in the experimental ^^Al samples positioned normally to 
the proton beam axis at the beam center and at distances 
of 40-430 mm from the beam axis. The 2'^Al(n,p)2^Mg 
cross sections were derived from the MENDL2 library 
|16[ . The distance-dependent <I>„/$p ratios measured via 
^Al(n,p)'^^Mg reaction are shown in Fig. ^ together with 
the results of the SSNTD method. 

The technique has also made it possible to detect 
2'^Mg together with ^^Na, ^^Na, and '^Be, when ^'^Al was 
irradiated in the proton beam. In this case, the neutron- 
to-proton flux density ratio in the beam is estimated as 



where (t_ 



N 



Be/^Na/'^Na 



N Mg 



'Be,^Na,Na 
,x 

77 "'Mg 



(20) 



'Mg 



a„ and a„ 



"Na 



are the neutron 



spectrum-weighted cross sections of the above reactions 
calculated by formula (|l^) using the excitation functions 
from MENDL2 database [0; aiy'',^^^'', and CTp^^ 
are the 27Al(p,x)24Na, 27Al(p,x)22Na, and 27Al(p,x)7Be 
reaction cross sections; n'*^", n''''^", n"*^^ , and N^-^^ 
are numbers of the ^^Na, '^'^Na,, ^^Mg, and ^Be nuclei 
produced in the Al samples with allowance for their 
decays under irradiation. 

The techniques described above were used in the 
experiments with proton beams of 0.07, 0.1, 0.13, 0.2, 
0.8, 1.0, 1.6, and 2.6 GeV. 

As seen from Figs. ^ and ^ the neutron component 
estimates obtained by the both techniques (SSNTD and 
7 spectrometry) are essentially the same and are small. 
It should be noted that the 2''8pb(p,x) reaction products 
measured at 1 GeV proton energy (see Table II) include 
none of the nuclides producible in the (n,xn) reactions 
up to X = 5 (2'''*™Pb), for such nuclides are either stable 
or long-lived. Our estimates have shown that, in the 
case of Pb isotopes of masses less than or equal to 
^''''Pb, additional production of daughter nuclei by (n,xn) 
reactions is below the level of experimental errors. 
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FIG. 2. The neutron-to-proton mean flux density ra- 
tios versus distance from the proton beam axis calculated 
using: (1) the SSNTD measurements (upper plot); (2) 
^'^Al(n,p)'^^Mg measurements (bottom plot). 

B. Monitor reactions 

Contrary to ||] , the present work uses the ^ ''Al(p,x)^^Na 
monitor reaction. Current practice suggests that three 
monitor reactions should be used on ^^Al, namely, 
^''Al(p,x)^**Na for short-term irradiations, ^'^Al(p,x)^^Na 
and ^^Al(p,x)^Be for long-term irradiations. 

Therefore, we have made additional experiments to 
study the ^*Na and ''Be yields relative to ^^Na. The 
^^Na cross section values were obtained from Since 
^^Na, ''Be, and ^^Na are produced in one and the same 
irradiated sample, their cross section ratios (represented 
by a ^"-/cr as an example) in a general form 

allowing for the neutron component can be calculated 
as 
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FIG. 3. The neutron-to-proton mean flux density ratios in 
the proton beam versus proton energy as calculated from the 
7-spectrometry data. 



'■NajBe 



A, 



^NajBe 



(rye) 



F 



^0 



Be 



(21) 



Since the <I>„/$j, ratio does not exceed ~ 2% at any 
proton beam energy (see Fig. ^), the following simplified 
formula can be used: 



/"NajBe 
22 N a 



A, 



^NajBe 



(rye) 



F' 



'Na 



^0 



(?ye) 



2'iNaJBe p^-^NaJBe 



(22) 



The urgent necessity for high-precision monitoring has 
made us to extend our measurements to the entire energy 
range used in our experiments, although a single energy 
I GeV) is treated in the present work. Table | and Fig. 
show the measurement results. 



C. 7-spectrometer detection efficiency allowing for 
the height factors 

The 7 spectra have to be measured for a broad range 
of heights H (the source-detector distances) so that 
early counting may begin without overwhelmingly large 
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TABLE I. The monitor reaction cross sections [mb] averaged over the experiments 



Proton energy, 


The ^^Al(p,x)^^Na 


The measured reaction cross sections. 


GeV 


monitor cross 


Shown in brackets 


are the errors disregarding/allowing 




sections used in 


for the ^^Al(p,x 


)^^Na reaction cross section error 




this work 


^'Al(p,x)"'*Na 


"'Al(p,x)'Be 


0.067 


A A 1 1/1 

24.4 ± 1.4 


11.3 ± ( 0.5 / 0.8 ) 


0.76 ± ( 0.20 / 0.21 ) 


0.097 


19.1 ± 1.3 


11.0 ± ( 0.3 / 0.8 ) 


0.97 ± ( 0.07 / 0.10 ) 


0.127 


17.0 ± 1.3 


10.1 ± ( 0.3 / 0.8 ) 


1.14 ± ( 0.06 / 0.11 ) 


0.147 


16.1 ± 1.2 


Q7_i_/'n/i / r\ ^ \ 


i/i/i_i_/'nii / 


0.197 


15.1 ± 0.9 


9.8 ± ( 0.4 / 0.7 ) 


1.48 ± ( 0.04 / 0.10 ) 


0.8 


15.5 ± 0.9 


12.7 ± ( 0.3 / 0.8 ) 


6.3 ± ( 0.3 / 0.4 ) 


1.0 


15.0 ± 0.9 


12.5 ± ( 0.8 / 1.1 ) 


7.5 ± ( 0.3 / 0.5 ) 


1.2 


14.6 ± 1.0 


12.8 ± ( 0.3 / 0.9 ) 


8.3 ± ( 0.2 / 0.6 ) 


1.4 


13.9 ± 1.0 


12.7 ± ( 0.4 / 1.0 ) 


8.9 ± ( 0.3 / 0.7 ) 


1.5 


13.5 ± 1.0 


12.8 ± ( 0.3 / 1.0 ) 


8.8 ± ( 0.3 / 0.7 ) 


1.6 


13.2 ± 1.0 


11.6 ± ( 0.3 / 0.9 ) 


8.9 ± ( 0.2 / 0.7 ) 


2.6 


11.7 ± 0.9 


10.6 ± ( 0.3 / 0.9 ) 


9.2 ± ( 0.2 / 0.7 ) 



counting rates, and later measurements do not suffer 
from too low rates. 

Therefore, the 7 spectrometer was efficiency calibrated 
in a 60 < E^< 2600 keV range for nine heights within a 
40-550 mm interval using the standard set of 7-radiation 
sources (^^Mn, 57c, ^OCo, ^Sy, i09cd, "^Sn, "Sga, 
"^Cs, i39Ce, i52eu, 228Th, 241 Am, and 22Na). The 
energy dependence of the detection efficiency obtained 
for a base height was approximated by the cubic spline 



SbaseiE) = exp 



■fc+3 



^P,.(lni?r 



(23) 



se ^base 

(E) X 



\ 



k+3 fc+3 

^^M,^i(lni?)''+^--2fc 

i—k j—k 



(24) 



where 



experimental efficiency of the spectrometer at the i-th 
energy E,. 

The radioactive reaction products of as long a half- 
life as possible were sought for and identified for each 
of the irradiated samples. Therefore, starting from a 
certain moment, the measurements were made at the 
lowest height (H=40 mm) , where the effect of the cascade 
summation of 7 quanta is negligible. The magnitude of 
this effect was determined experimentally. 

An optimal schedule for changing the heights to 
preserve an approximate constancy of the maximum 
possible spectrometer load was very difficult to select 
beginning from the very start moment of measuring a 
sample irradiated. Therefore, the height dependence of 
the detection efficiency was analyzed for each of the 7 
energies of the calibration sources. The analysis has given 
the dependence expressed as 



eiE,H) 



'-hast 



{qi +q2 - ItlE + Hi 



base 



(gi +92 - InE + H) 



(26) 



1, for ln£' < InE'o, = 300 keV, 

5, for ln£' > InE'o, 



J1+J2 

i=l 



'base 



(Ei) - Ebase (Ei) 



(25) 



where Sbase{Ei) and Ael^^^{Ei) are the value and its 
absolute error for the calculated efficiency of the spec- 
trometer at the i-th energy Ei; Pi are polynomial coeffi- 
cients; M~j^ is the covariance matrix of the polynomial 
coefficients; L is the number of degrees of freedom 
when determining x^- _|_ jg the total number of 

experimental points; Sbasei-^i) ^^'^ '^^tasei^i) 
experimental values and their absolute errors for the 



Ae{E,H)=e{E,H) 



Ae, 



base 



£ba 



-A{H-Hbase?-X 



{^^EY)-xl 



[qi +q2-\nE + Hbasef -[qi+qi-lnE + Hf 



(27) 



where e {E, H) and Ae{E, H) are, respectively, the value 
and the absolute error for the dependence of the spec- 
trometer detection efficiency on energy and height; qi 
and 52 are the fitting parameters; V~^^ are the covariance 
matrix elements of the parameters qi and 92; X5 is given 
by Eq. (H). 
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FIG. 4. The ■^^Al(p,x)^^Na monitor cross sections mea- 
sured in this work and from previous works: MI85 - [18], 
MI86 - [19], MI89 - [20], MI90 - [21], MI93 - [22] MI95 - [7] 
MI96 - [23], MI97 - [24]. 



From formula ( p6[ ) it is readily seen that any of the 
above mentioned heights can be selected, but the extreme 
heights are not desirable, for they involve a rise of 
the systematic error when renormalizing the rest of the 
heights to the basic height. If, however, the height of the 
prevailing measurements is selected, the systematic error 
decreases during the renormalization. 

Fig. H shows the detection efficiencies measured at 
different heights (the dots). The curves in the figure 
present the dependence of the efhciencies on height and 
energy as inferred from formula (p6|). 

The analysis based on formulas (^) and ( |27| ) has 
shown that the systematic errors due to uncertainties in 
the detection efficiency do not exceed 5%. 
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FIG. 5. The measured and simulated height-energy detec- 
tion efficiency of the CANBERRA GC2518 Ge-detector-based 
7 spectrometer for different heights of the 7 source. 



for permitting the results of tentative processing to be 
analyzed again after the computer-aided consolidated 
processing of the 7 spectra by interactive fitting of the 
7 peaks. Namely, we find out whether the peak regions 
are either multiplets, real peaks that fail to satisfy the 
search criteria, spurious peaks, etc. 

The interactive processing has made it possible to 
improve the accuracy and reliability of analysis of the 
measured 7 spectra, especially in the case of poorly 
resolved spectra with poor statistics. 

The processed 7 spectra are then united to form the 
input file for the ITEP-developed SIGMA code, which 
plots the time variations of a selected 7-line intensity 
and makes use of the energy and the half-life of that 7 
line to identify the product nuclide using the PCNUDAT 
database |Q and the Firestone's Table of Isotopes H^] 
and to determine the product nuclide cross sections using 
formulas 



D. Processing of 7 spectra and cross section 
determination 

The measured 7 spectra of irradiated ^°®Pb are ex- 
tremely complicated due to the numerous 7 lines. De- 
spite using the Ge detector with its ultimate attained 
resolution, the spectra still include many multiplets, so 
the GENIE-2000 code [^ which is designed to process 
complicated 7 spectra has been used. The code is notable 



IV. MEASURED PRODUCT YIELDS FROM THE 
208p3(p X) REACTION 

Table II presents the results of determining the reac- 
tion product yields in 1 GeV proton-irradiated ^°^Pb. 
Out of 114 yields presented, 8 are independent yields of 
ground states (i), 15 are independent yields of metastable 
states (i(I]TOj)), 15 are independent yields of metastable 
and ground states (i(STOj-l-g)), 65 are cumulative yields 
(c), and 11 are supra cumulative yields, when the 
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correction may exceed the determination error (c*). It 
should be noted that in case the (c) values are to be 
reconstructed from (c*) using of formula (15), we must 
have the crj"'" values that have yet to be determined, but 
are expected to be measured in the future. As mentioned 
above, however, the simulation codes can be tested using 
the (c*) values. 

^From Table II we see that the experimental errors 
range from 7% to 30%. The tabulated values are 
calculated as follows. Since most of the results are 
obtained by averaging over a set of ((7^ ± Act.;) values 
calculated from different 7 lines, the mean and the 
experimental errors were calculated as 

a^^^, where W,^l/^aj. (28) 

The A(7i were determined using the error propagation 
formulas 1231 



Aa 



(n-l)E.W, 



Act = 



(29) 



(30) 




In the calculations, the largest of Aa and Act was 
taken to be the experimental error Act. Allowing for 
the monitor error, the total error in the measured cross 
sections was calculated as 



(31) 



Our analysis has shown that the main contribution to the 
total error is from uncertainties in the nuclear data (the 
absolute gamma yields and cross sections of the monitor 
reactions) . 



V. COMPARISON WITH EXPERIMENTAL 
DATA OBTAINED ELSEWHERE 

In addition to the nuclide production cross sections 
measured in this work. Table II also presents data from 
other works, namely, the results obtained by the "inverse 
kinematics" technique (^"^^Pb bombarding at GSI ||]) 
and the results of measuring the nuclide production in 
1 GeV proton-irradiated "°*Pb obtained by a technique 
that is similar to ours by the group of R. Michel at 
Hannover University, Zcntrum fiir Strahlenschutz und 
Radiookologie (ZSR) g. 

It should be noted that the GSI data are measured only 
as independent yields, so they are summed up for their 
isobaric chains to be compared with the respective cumu- 
lative yields determined at ITEP and ZSR. We estimated 



errors for the GSI cumulative yields as square roots of 
sums of squared errors of all independent yields con- 
tributing to the corresponding cumulative cross sections 
and also took into account their systematic errors which 
vary from 9 to 25% . The GSI technique does not dis- 
tinguish if a product nuclide is in a ground or metastable 
state, hence, we do not compare our data with the GSI set 
in the case we measured only metastable or ground states 

^204mp]-| 197inpj_| 198mlrpj 196mrpj 197mjjg 195mjjg 

i93mjjg/i98mAu^ 1989 Au, i^Sffir, is^^Os, ^^^""Re, ^^^""Te, 

119mr£Q 120mg^ 114mj^ llOm^g 106m^g lOlniJ^J^ dOm-y 

and ^^"^Rh) as well as when there is a transition of a 
metastable state to a product out of the given decay chain 
(198T1, i90lr, i52xb, i49Gd, ^^iXe, 96tc, 95Tc, 95Nb, 
s^Nb, s^Zr, s'^Y, and ^^Sr). 

Let us mention here that the ZSR work [E9[ was carried 
out using natural lead, rather than high Pb-enriched 
lead used at ITEP and GSI. We think, it is this fact 
that can be primarily responsible for some divergences 
between the ZSR data and the GSI and ITEP ones. The 
importance of this circumstance is indirectly confirmed 
by a 1.5 times higher mean-squared deviation of the ZSR 
data from the GSI data compared with the mean-squared 
deviation of the ITEP data from the same GSI data. The 
mean-squared deviations between two data sets (1 and 
2) calculated as < Fi_2 >= 100% x (10^ - 1) (where 
D =< log{au/a2if > ) proved to be 48% and 32% for 
ZSR-GSI and ITEP-GSI comparisons, respectively. 



VI. SIMULATION OF EXPERIMENTAL DATA 

The present work is aimed at determining the nuclear 
cross sections to be used in designing ADS facilities. An 
attempt to obtain the necessary nuclear cross sections 
only from experiments would involve impractical levels 
of expense and effort. Therefore, simulation techniques 
must be used for that purpose. The simulation approach 
has the advantage that it also can be used for many other 
purposes. On the other hand, the current accuracy and 
reliability of simulated data are inferior to experiment. 
Also existing simulation codes have different predictive 
abilities when used to study the reactions that are of 
practical importance. 

Considering this, the present work is aimed also at 
verifying the simulation codes used most extensively for 
this purpose in order to estimate their predictive abilities 
and to stimulate efforts to improve them. 

The following seven simulation codes were examined 
to meet these requirements: 

• the CEM95 cascade-exciton model code |3Ct| , 

• the latest version of the improved cascade-exciton 
model code, CEM2k, 

• the CASCADE cascade-evaporation-fission trans- 
port code [|3|, 
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TABLE II. Experimental product nuclide yields [mb] in 1 GeV proton-irradiated '^"^Pb compared with both recent GSI data 
measured in inverse kinematics [8] and the ZSR measurements on ""*Pb [29] 



Product 


Ti/2 0, 

1 ±Lz — H \L- — 1 U 


Type 


Yield (this work) 


Work il 
1 a_t — 1 


Work S 

1 , 1 V 




6.243d 




4 60 + n 29 






205 g; 


15. Sid 




6 20 ± 40 




7 09 ± 90 


204 g J 


11.22h 


ifml4-m2+el 


5.29 ± 0.80 




6.03 ± 0.95 


203 g; 


11.76h 




4 84 ± 59 






204mp^ 


67.2m 


i(m) 


11.0 ± 1.0 






203 


51 873h 


Q 


31.5 ± 2.1 


28.7 ± 3.1 




201 


9.33h 


c* 


26 9 + 2 4 


20.4 ± 1.9 




200 p^ 


21. 5h 


Q 


18.2 ± 1.2 


18.2 ± 2.0 


27 8 ± 3 5 


198 p^ 


2.4h 




8.9 ± 2.1 


14.0 ± 1.3 




197m p^ 


43m 




17.9 ± 4.0 






202rpj 


12.23d 


Q 


18.8 ± 1.2 


40.0 ± 4.0 


22.0 ± 2.7 


201 rpj 


72.912h 


Q 


43.7 ± 2.9 


37.3 ± 3.7 


53 5 ± 6 6 


200rpj 


26. Ih 


c 


40.6 ± 2.6 


35.2 ± 3.7 




200rpj 


26. Ih 


i ( TV| -t-P' 1 


22.7 ± 1.5 


17.0 ± 1.7 


22.3 ± 6.1 


igg-pj 


7.42h 


c 


38.5 ± 5.2 


34.3 ± 3.4 




198ml rpj 


1.87h 


i(mH-m2) 


17.6 ± 3.6 






198rpj 


5.3h 


Q 


35 9 ± 5 






196m rp^ 


1.41h 


ifm) 


34.8 ± 4.4 






203He 


46.612d 


Q 


4 03 ± 27 




3 66 ± 45 




23. 8h 


ifml 


10.7 ± 0.7 






195mTT 


41. 6h 


ifml 


13.6 ± 2.0 




13.3 ± 1.8 


193mTT 


11. 8h 


i(m) 


18.9 ± 2.5 




10.8 ± 2.3 




4.85h 


c 


35.2 ± 2.8 


31.3 ± 3.4 




198m^^ 


2.27d 


ifm) 


1.01 ± 0.14 




1.25 ± 1.11 




2.69517d 


]( m-l-ffl 


2.11 ± 0.22 


1.96 ± 0.23 




"*Au 


2 695 17d 




1 09 ± 30 








6.183d 


ifml4-m2+el 


4.13 ± 0.35 


4.02 ± 0.47 


3.88 ± 0.47 


An 


186.098d 


c 


48.7 ± 5.5 


28.4 ± 3.3 


51.1 ± 6.6 


"■•Au 


38.02h 


i('ml-l~m2+e"l 


7 06 ± 75 


6 33 ± 75 


6 85 ± 92 


"2 An 


4.94h 




46 9 ± 6 6 


39 9 ± 4 6 






4.94h 


ii'ml -^m24'e'l 


11.6 ± 1.7 


9.2 ± 1.1 




191p^ 


2.802d 


Q 


41.8 ± 4.2 


44.4 ± 5.5 


37.9 ± 4.8 


189 p^ 


10.87h 


Q 


46.8 ± 4.8 


40.4 ± 5.0 




188 p^ 


10. 2d 


Q 


40.5 ± 2.9 


38.4 ± 4.7 


42.8 ± 5.4 


186 p^ 


2.08h 


c* 


33 5 ± 2 3 


32.9 ± 4.1 




190 


11.78d 




69 ± 06 






188 jj. 


41. 5h 


Q 


43.2 ± 3.2 


40.9 ± 5.4 




188 jj. 


41. 5h 




2 93 ± 69 


2 48 ± 33 




186 jj. 


16.64h 


\ 


20.8 ± 1.9 




22.5 ± 3.1 


185 jj. 


14. 4h 


c* 


34.8 ± 2.3 


39.4 ± 5.2 


39.4 ± 7.9 


184 


3.09h 


c* 


39 5 ± 3 


36 9 ± 4 8 






93. 6d 


Q 


41.8 ± 2.8 


38 1 ± 5 3 


43 ± 5 3 


183m Qg 


9.9h 


c 


23.2 ± 1.5 








22.10h 


Q 


42.0 ± 2.8 


34.2 ± 4.8 




Re 


70. Od 


Q 


41.7 ± 2.9 


36 3 ± 5 3 


38.2 ± 4.8 


182m pg 


12. 7h 


Q 


45.2 ± 3.7 






i^^Re 


19. 9h 


Q 


43.1 ± 5.9 


37.0 ± 5.4 


45 9 ± 5 9 


i™Re 


19.5m 


c* 


48.2 ± 4.2 


44.7 ± 6.6 






135m 


C 


30.1 ± 3.5 


23.4 ± 3.6 






2.5h 


c 


28.0 ± 3.9 


29.0 ± 4.5 




i^^Ta 


8.09h 


c 


35.0 ± 3.6 


28.8 ± 4.7 




i^^Ta 


3.14h 


c 


30.9 ± 3.9 


26.3 ± 4.3 




"^Ta 


36.8m 


c* 


17.3 ± 2.3 


27.4 ± 4.5 




175 jjf 


70d 


c 


31.3 ± 2.3 


28.3 ± 4.8 


34.1 ± 4.1 


173 jjf 


23.6h 


c 


28.4 ± 2.6 


25.2 ± 4.3 


39.0 ± 4.9 


172 jjf 


1.87y 


c 


24.1 ± 1.6 


24.6 ± 4.2 


24.4 ± 3.1 
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Product 


Ti/2 


Type 


Yield (this 


work) 


Work [ 






Work 1 


29 






12. Ih 


c 




99 Q 4_ o q 
zz.y ziz o.y 










170 jjf 


le.Olh 


c 


99 1 -4- fi B 

ZZ. L U.O 


90 q 4_ o c 
zu.o ZIZ o.o 


21.2 zb 3.0 




6.70d 


c 


9'3 q -1- 1 7 


9A 7 -\- A A 

Z'^:. ( ZIZ 










Lu 


6.70d 


i(ml+m2+g) 


n 1 Q + n 


n 1 oq 4_ f) nqy 

U.-LOO ZIZ U.UOi 












8.24d 


c 


9fi 1 -1- 1 S 


1 f:; 4_ on 
ID.D zt O.U 


31.3 zb 3.9 


"°Lu 


2.012d 


c 


21.7 ± 2.9 


20.9 zb 3.7 












34.06h 


c 


18.6 ± 1.2 


12.1 zb 2.2 


26.4 ± 3.7 


169Yb 


32.026d 


c 


20.9 ± 1.5 


18.1 zb 3.4 


24.3 zb 3.0 


166Yb 


56.7h 


c 


16.1 ± 1.1 


13.7 zb 2.6 


16.4 ± 2.3 




9.25d 


c 


19.4 ± 4.0 


14.0 zb 2.7 


21.2 ± 2.6 




30.06h 


c 


14.4 ± 1.4 


13.3 zb 2.6 












28.58h 


c 


o.o 


_1_ 

Hz 


0.6 


7 9 
i .Z 


4- 
zt 


1 

















8.14h 


c 


O. 1 o 


-L 
ZC 


0.45 


O.U 


4- 

ZIZ 


1 




1 
1 












9.9h 


c* 


O.DD 


It 


u 


97 
Z / 


9 Sfi 
Z.oD 


4_ 
zt 


0. 


DO 












5.32d 


c 


/I 1 


-4- 


U 


oy 


9 79 
Z. 1 z 


4- 
zt 


0. 


fi9 
DZ 


K f^9 


zb 0.70 


153 


2.34d 


c* 


Z.OZ 


_l_ 
zt 


U 


9^; 
zo 


9 /10 


4_ 
zt 


0. 


KA 
04 


9 ^^ 
Z.Ol 


zb 0.40 


152 rp^ 


17.5h 


c* 


1 n 


_l_ 
zt 


U 


1 7 

1 ( 






















240.4d 


c 


Z.OD 


_l_ 
zt 


U 


9/1 
z4 


9 1 S 
Z. lo 


4_ 
zt 


0. 


01 


Q 1 n 
o.iU 


zb 0.38 


i«Gd 


9.28d 


c 


9 9/1 


It 


U 


lo 












o.UO 


zb 0.38 


"«Gd 


48.27d 


c 


1 9^^ 
l.ZD 


zt 


U 


HQ 

uy 


1 99 
l.ZO 


4_ 
zt 


0. 


90 

zy 


l.Oo 


zb 0.21 




24. Id 


c 


u.yo 


_|_ 

zt 


U 


oU 


IIS 
1. lo 


4_ 
zt 


0. 


9Q 

zy 


i.y / 


zb 0.29 


i«Eu 


4.61d 


c 


1 fi9 


_|_ 
Zt 


U 


1 9 
IZ 


1 1 7 
1. 1 i 


4- 
zt 


0. 


9S 
Zo 










"«Eu 


4.61d 


i 


U.O / 


_l_ 
zt 


U 


uo 


O 1 SI 
U. lol 


4_ 
zt 


0. 


0/17 
U4 / 










i«Pm 


265d 


c 


1 09 


_l_ 
zt 


U 


lo 


o sp; 

U.oO 


4_ 
zt 


0. 


99 
zZ 


1 on 
l.UU 


zb 0.13 


«9Ce 


137.640d 


c 


yj.oo 


4- 
it 


U 


UD 


n /I p; 


4- 
zt 


0. 


10 


n C9 
U.oz 


± 0.10 


121m 


154d 


i(m) 


n AA 


_L 
Zt 


U 


U-l 












U.OO 


zb 0.07 


Te 


19.16d 


c 


l.U / 


It 


U 


1 1 












U. / y 


zb 0.10 


119mrp^ 


4.70d 


i(m) 


n AH 


It 


U 


0/1 
U4 




















120msb 


5.76d 


i(m) 


U.04 


_|_ 
zt 


u 


op; 
Uo 












U.OO 


± 0.07 


"-""In 


49.51d 


i(ml+m2) 




u.yo 


_l_ 
zt 


u 


ly 












l.U < 


zb 0.16 


llOm^g 


249.76d 


i(m) 


111 
1.11 


zt 


U 


OQ 

uy 












1 Q9 
l.OZ 


zb 0.17 


106m^g 


8.28d 


i(m) 


fiQ 


It 


u 


OS 

Uo 












o no 

u.y/ 


zb 0.14 




41.29d 


c 


U.DD 


It 


u 


1 9 
Iz 


O 7/1 
U. / 4 


4_ 
zt 


0. 


17 


i.U4 


± 0.14 


Rh 


35.36h 


c 


4. Do 


4_ 
zt 


U 


KA 
04 


o. 10 


4_ 
zt 


0. 


51 










i°i™Rh 


4.34d 


c 


1 9Q 

i.zy 


4_ 
zt 


U 


ID 




















103 Ru 


39.26d 


c 


Q ft/1 


4_ 
zt 


u 


9fi 
ZD 


Q OQ 
o.UO 


4_ 
zt 


0. 


50 


/111 
4.11 


zb 0.53 


96Tc 


4.28d 


i(m+g) 


1 90 
1 .zu 


4- 
zt 




u 


OQ 

uy 












1 /10 

i.4y 


zb 0.19 


95^^ 


20.0h 


c 


l.Oo 


4_ 
zt 


U 


lo 




















s^Nb 


23.35h 


i 


9 "^1 
Z.Ol 


4- 
zt 





19 


Z. lO 


4- 
zt 


0. 


34 










^^Nb 


34.975d 


c 




4- 
zt 





34 




















95 Nb 


34.975d 


i(m+g) 


3.03 


zt 





20 












o.Oo 


zb 0.56 


«5Zr 


64.02d 


c 


2.34 


zb 





15 


1.58 


zb 


0. 


28 


2.32 


zb 0.29 




78.41h 


c 


2.30 


zb 





16 












2.82 


zb 0.35 


8«Zr 


83.4d 


c 


0.76 


zb 





08 


0.97 


zb 


0. 


15 


1.19 


zb 0.15 


90m Y 


3.19h 


i(m) 


4.82 


zb 





39 




















88y 


106.65d 


c 


4.03 


zb 





27 


3.72 


zb 


0. 


58 










88y 


106.65d 


i(m+g) 


3.41 


zb 





25 


2.76 


zb 


0. 


44 


3.74 


zb 0.46 
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c* 
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zb 
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»5Sr 
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c 
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zb 
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»«Rb 
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± 





66 
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± 
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c 
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zb 





28 


2.82 


± 


0.45 


3.96 


zb 0.49 


82m 


6.472h 


i(m) 


2.73 


zb 
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35.30h 


i(m+g) 


2.17 


zb 





14 


1.55 


zb 


0.24 


2.62 


zb 0.50 




119.779d 


c 


1.33 


zb 





09 


1.18 


zb 


0.19 


1.61 


zb 0.20 




17.77d 


i 


1.86 


± 





18 


1.66 


± 


0.27 


2.24 


± 0.28 


59pg 


44.472d 


c 


0.91 


± 





08 


0.69 


± 


0.11 


1.05 


± 0.14 


«5Zn 


244.26d 


c 


0.79 


zb 





19 


0.42 


zb 


0.07 


0.66 


± 0.17 




83.79d 


i(m+g) 


0.35 


zb 





06 












0.37 


± 0.05 
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• the INUCL 
fission code [t 



cascade-preequilibrium-evaporation- 



• the LAHET (both ISABEL and Bertini options) 
cascade-preequihbrium-evaporation- fission 
transport code [ p5[ , 

• the YIELDX semi-phenomenological code , 

• the CASCADE/INPE cascade-evaporation- 
preequihbrium- fission-transport code \ ^7\ . 

All the codes, except for CASCADE/INPE and 
CEM2k are described in some detail in |^ . 

The intranuclear part of the CASCADE/INPE code 
package based on the Dubna cascade model |3|], is 
used to simulate the characteristics of projectile interac- 
tions with target nuclei. Recently the code was upgraded 
at the Institute of Nuclear Power Engineering (Obninsk, 
Russia) Principal modifications of the code include: 

• a special routine for calculating the precompound 
spectra of a particles, which have been demon- 
strated to play an important role in the produc- 
tion of long-lived radioactivity in different heavy 
targets; 

• the code was upgraded to allow for the description 
of angular distributions of nucleons based on the 
Kalbach parametrization for projectile energies be- 
low 0.8 GeV; 

• a new realization of the Weisskopf evaporation 
approach used in |Q has been introduced. This is 
based on inverse-reaction cross-section calculations 
using various optical potentials appropriate for 
various mass and energy regions; nuclear level 
densities are calculated taking into consideration 
nucleon pairing to evaluate the excitation energy 
of residual nuclei; 

• the latest version of the mass table of the nuclides 

has been used for binding-energy calculations. 

The improved Cascade-Exciton Model (CEM) of nu- 
clear reactions was developed and incorporated in the 
code CEM97 in the Theoretical Division of Los Alamos 
National Laboratory, as an improvement to the code 
CEM95 Q. It is described in detail in therefore 
we will not elaborate here. CEM2k is a next step in the 
improvement of the CEM; it differs from CEM97 mainly 
in the details of the transitions from the cascade stage 
of a reaction to the preequilibrium one, and from the 
latter to equilibrium decay. This preliminary version of 
CEM2k has less preequilibrium emission than the earlier 
versions. The changes were motivated by discrepancies 
with the recent GSI data ^ in the earlier versions of 
the model. It is briefly surveyed in [|2|, it is still under 
development and will be described in a future paper. 

Contrary to the simulated data, the experimental 
results include not only independent, but also (and 



mainly) cumulative and supra cumulative, residual prod- 
uct nuclei. To get a correct comparison between the 
experimental and simulation data, the cumulative yields 
must be calculated on the basis of the simulated inde- 
pendent yields. If the production chain of n radioactive 
nuclei is presented as 



i 
1 



i 
2 



V2 



(33) 



(where v\,...,Vn-\ are the branching ratios of the 
respective nuclides), the simulated cumulative and supra 
cumulative yields of the n-th nuclide can be calculated 

as 



n-l 

E 

i=\ 



n 



(33) 



cum ind 

CT„ — cr„ 



And 



ind 



i=l 



n 

0=^ 



(34) 



The branching ratios of the decay chains were taken 
from 1^ , considering that the branched (due to isomeric 
transitions and a decay) isobaric chains can always be 
presented as a superposition of linear chains. 

To get a correct comparison between results obtained 
by different codes, the calculations were renormalized 
to the same cross sections for proton-nucleus inelastic 
interactions. We calculated these cross sections by the 
semi-empirical formula |40[| 



a.nei = ^hA'>-^f{A)g{E), ( in mb) 

f{A) = 1 + 0.016 sin(5. 3-2.63 log ^), 

g{E) = 1 - 0.62exp(-£;/200) sin(10.9i; - 0.28), 



(35) 



where A is the mass number of the target and E is the 
energy (in MeV) of the projectile proton. 

The cross section for the p-l-^^'^Pb inelastic interactions 
as calculated by formula ( [Sq ) is of 1857.8 mb at 995 MeV, 
the incident energy of our experiment. 

If an experiment-simulation difference of no more 
than 30% (0.77< (Tcaic/<^exv <l-3) is taken to be the 
coincidence criterion Bl| , the simulation accuracy can be 
defined as the ratio of the number of such coincidences 
to the number of the comparison events. The 30% level 
meets the accuracy requirements of cross sections for 
nuclide production to be used in designing ADS plants, 
according to The mean simulated-to-experimental 

data ratio can be used as another coincidence criterion 



< F > = 10 



exp,t 



(36) 
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with its standard deviation 

S'(< F >) = < (log [(JcaLi/(Texp,i) - log(< F >)f >, 

(37) 

where <> designates averaging over all the experimental 
and simulated results used in the comparisons {i = 
l,...,Ns). 

The mean ratio < F > together with its standard 
deviation S'(< F >) defines the interval [< F > /S{< 
F >) , < F > xS{< F >)] that covers about 2/3 of 
the simulation-to-experiment ratios. A logarithmic scale 
is preferable when determining the factor < F > rather 
than a linear scale, because the simulation-experiment 
differences may be as high as a few orders of magnitude. 

We apply the above two criteria together with our 
results shown in Figs. |^-||, to infer conclusions about the 
predictive power of a given code. 

The default options were used in all of the simula- 
tion codes without modifying the codes to get optimal 
agreement with the data. All the calculations were 
made before any experimental results were obtained, 
except the results from CEM2k. With such an approach, 
our comparisons demonstrate the real predictive power, 
rather than the descriptive power of the codes. 



VII. COMPARISON OF EXPERIMENT WITH 
SIMULATIONS 

The results obtained with the codes are presented in 

• Fig. ^, which shows the results of a detailed 
comparison between simulated and experimental 
independent and cumulative products (criterion 1); 

• Fig. 0, which shows the results of a detailed 
comparison between simulated and experimental 
independent products of all isotopes of Tm, Ir, 
and Tl measured in this experiment (black squares) 
together with the data obtained by the reverse- 
kinematics method at GSI (black stars) 

• Fig. ^, which shows the statistics of the simulated- 
to-experimental data ratios (criterion 2); 

• Fig. 1^, which shows the simulated mass distribu- 
tions of the products together with the measured 
cumulative and supra cumulative yields of nuclides 
that are in immediate proximity to the stable 
isotope of a given mass (the sum of such yields from 
either side in cases when both left- and right-hand 
branches of the chain are present). Obviously, the 
simulation results do not contradict the experimen- 
tal data if calculated values run above the experi- 
mental data and follow a general trend of the latter. 
This is because direct 7 spectrometry identifies only 
the radioactive products, which generally form a 
significant fraction of the total mass yield but are 



never equal to the total mass yield when a stable 
isobar is produced. 



193mjjg 
110m ; 



Table III presents quantitative information concerning 
the agreement of the simulated yields with experimental 
data for each of the simulation codes, namely: 

• the total number of measured yields, Nj^; 

• of them, the number of the measured yields selected 
to compare with calculations, Ng. We reject the 
following nuclides from our comparison in the cases 
where: 

1. The measured product is of metastable or 
just ground state, namely, ^°'''"Pb, ^^'^'"Pb, 

198mlrp2 196mrp^ 197mjjg 195mjjg 
198m Au^' 198ffAu, ^""^Hl,' 183™Os, 

i2i™Te, "9"Te, ^'^^"'Sh, "'^"Tn, ""™Ag, 
106m loi^Rh, 90™Y, and 82™Rb; 

2. There is a transition of a metastable state to a 
product out of the given decay chain, namely, 
198T1, i90lr, 152Tb, i49Gd, i2iTe, ^^Tc, ^^Tc, 
95Nb, ^5^h, s^Zr, s^Y, and ^^Sr; 

3. There is a strong correlation between a mea- 
sured cumulative yield decaying into another 
one, namely, "^^a — > i^sjjg^ ivonf — , 
i™Lu, i69Lu i69Yb, and ^^STb ~^ ^^^Gd. 
The cumulative yields of the precursors in 
all the above chains are almost equal to the 
cumulative yields of the daughters, that is 
why only the daughter yields were used in 
our comparison, to prevent double counting. 
Also, in case of a strong correlation between 
the cumulative and independent yields of a 
product (®*Y), only the independent yield was 
used for comparison. 

• of them, the number of the product nuclei whose 
yields were simulated by a particular code, Ng; 

• the number of comparison events when the sim- 
ulated results differ from the experimental data 
by not more than 30%, 31 a-nd the number of 
comparison events when the calculations differ from 
data by not more than a factor of 2.0, N^^ g-, 

• the mean squared deviation of the simulated results 
from experimental data, < F >, and its standard 
deviation, S{< F >). 

Since about a third of all secondary nuclei from our 
reaction are not spallation products, the ability of codes 
to simulate high-energy fission processes is an important 
criterion for their ability to work when the target is heavy 
enough to fission. Among the codes used here, LAHET, 
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FIG. 6. Detailed comparison between experimental and simulated yields of radioactive reaction products. The cumulative 
yields are labeled with a "c" when the respective independent yields are also shown. 
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FIG. 7. Isotopic mass distribution for independent products of Tm, Ir, and Tl isotopes. Black squares are our measurements, 
while filled stars show GSI data obtained in reverse kinematics. Results from different codes are marked as indicated. 
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FIG. 8. Statistics of the simulation-toexperiment ratios (criterion 2). 
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FIG. 9. The simulated mass distributions of reaction products together with the measured cumulative and supra-cumulative 
yields. Black line shows GSI data in reverse kinematics. 
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TABLE III. Statistics of comparison between experimental 
and simulated yields in 1.0 GeV proton-irradiated ^"^Pb 



Code 


1 Nt = 114, Ng = 76 1 


Nc,3/ 
Nc..o/Ns 


<F> S{<F>) 


LAHET-ISABEL 


39/59/76 


1.92 1.73 


LAHET-Bertini 


33/54/76 


2.04 1.71 


CEM95 


30/46/55 


2.16 2.00 


CASCADE 


29/55/72 


2.18 1.88 


CASCADE/INPE 


32/56/70 


1.85 1.59 


INUCL 


24/40/73 


2.89 2.15 


YIELDX 


27/49/76 


2.76 2.24 


CEM2k 


33/54/60 


1.62 1.45 



CASCADE, INUCL, CASCADE/INPE, and YIELDX 
simulate both spallation and fission products. The 
CEM95 and CEM2k codes simulate spallation only and 
do not calculate the process of fission, and do not provide 
fission fragments and a further possible evaporation of 
particles from them. When, during a Monte Carlo 
simulation of the compound stage of a reaction these 
codes encounter a fission, they simply remember this 
event (that permits them to calculate fission cross section 
and fissility) and finish the calculation of this event 
without a subsequent calculation of fission fragments. 
Therefore, results from CEM95 and CEM2k shown here 
reflect the contribution to the total yields of the nuclides 
only from deep spallation processes of successive emission 
of particles from the target, but do not contain fission 
products. This is explicitly reflected in a smaller number 
of the products simulated (the quantity Ns in Table III 
and in the shapes of the simulation curves in Figs. ^ 
and ^). To be able to describe nuclide production in 
the fission region, these codes have to be extended by 
incorporating a model of high energy fission (e.g., in the 
transport code MCNPX |4|], where CEM97 is used, it is 
supplemented by the RAL fission model [^). 

The following conclusions follow from the analysis of 
the results presented in Table III and in Figs. 6-9: 

1. All codes can reasonably adequately simulate the 
weak spallation reactions (the A>180 products), 
with simulation results differing from experimental 
data, usually within a factor of 2. The excep- 
tions are the yields of ^''•^Hg (underestimated by 
an order of magnitude by YIELDX) and i^Sjj. 
(underestimated by an order of magnitude by 
CASCADE/INPE). It should be noted that all the 
codes simulate a number of products (^^''ir, ^^^Tl, 
199T1, 200Tl(cum), 201T1, and ^ospb), which differ 
from experiment by less than 30%. 

2. In the range of the deep spallation reactions 
(150<A<180), the simulation codes have very dif- 



ferent predictive powers, namely, 

• the LAHET, CASCADE/INPE, and YIELDX 
predictions are very close to the experimen- 
tal data; the only exception is ^^^Lu(ind), 
whose measured yield is about a factor of 
two smaller than LAHET and YIELDX give 
and about a factor of two higher than the 
CASCADE/INPE prediction; 

• the CASCADE code simulates the A>160 
product yields quite adequately, except for 
-'^^^Lu(ind) whose measured yield is about ten 
times above the simulated value; however, as 
the atomic number of the product decreases 
below 160, we observe an underestimation of 
data by CASCADE; the deviation tends to 
increase with decreasing A (up to a factor of 

5); 

• the INUCL code underestimates these reac- 
tion products systematically by a factor of 
2-10 with the discrepancy increasing with 
decreasing A; 

• the new code CEM2k was developed taking 
into account the recent GSI measurements 

In the spallation region it agrees best 
with the data compared to all the other 
codes, although like its predecessor CEM95, it 
does not contain explicit treatment of fission 
fragments and should be supplemented in a 
transport code by a model of fission-fragment 
formation, or should be improved further by 
the development and incorporation of its own 
model of fission-fragment formation. 

3. In the range of fission products (40<A<150), the 
INUCL code predictions are in the best agreement 
with the data. As a rule, the INUCL-simulated 
yields differ from measured data by a factor of 
less than 1.5, except for ^^Fe, ^'^^Ce, and ^^'^Pm. 
The LAHET-simulated yields underestimate data 
by a factor of 1.5-10, except for ^®Zr and ^°^Ag. 
The specific predictions of isotope production cross 
sections (Fig. 6) of the semi-phenomenological 
code YIELDX both under- and overestimate the 
fission product data by a factor of up to 30, 
without showing any obvious patterns in the dis- 
agreement. In contrast, the YIELDX isobar cross 
sections are all greatly ovcrpredictcd, as shown 
in Fig. |. The CASCADE/lNPE-simulated yields 
of the 130<A<150 products are strongly under- 
estimated (up to 1-2 orders of magnitude), while 
the rest of the simulated fission product yields 
agree with the experimental data generally within 
a factor of 2. As a rule, the agreement of all 
codes with the data in the fission product region 
is worse than in the spallation region; therefore, 
development of a better model for fission-fragment 
formation is welcomed for any code. 
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VIII. CONCLUSION 

The interest shown m both the possible transmutation 
of nuclear wastes and the Spallation Neutron Source 
(SNS) facilities encourage us to anticipate that the 
accumulation and analysis of nuclear data for ADS and 
SNS applications will have the same growth in academic 
interest and practical commitments as was the case 
for nuclear reactor data during the last five decades. 
Therefore, experimental data on the yields of proton- 
induced reaction products as applied to the ADS and 
SNS main target and structure materials are of great 
interest and importance. It should be emphasized that 
the charge distributions in the isobaric decay chains 
are important as well. The information thus obtained 
would make it possible, first, to raise the information 
content of the comparisons between experimental and 
simulated data and, second, to reduce the uncertainties 
in experimental determination of the cumulative yields 

by establishing unambiguous relations between cr"^"™ and 
* 

^cum £qj, uiQj^y of the reaction products. 

We have measured in the present work 114 cross 
sections of nuclides produced in interactions of 1 GeV 
protons with ^"^Pb, of which 8 are independent yields of 
ground states, 15 are independent yields of metastable 
states, 15 are independent yields of metastable and 
ground states, 65 are cumulative yields, and 11 are supra 
cumulative yields. We have compared our data with 
previous measurements and with predictions of seven 
different codes used in many current applications to 
understand qualitatively and to estimate quantitatively 
their predictive powers. 

Regarding the codes benchmarked here, we conclude 
that none of them agree well with the data in the 
whole mass region of product nuclides and all should be 
improved further. In addition, the predictive power of all 
codes for data in the fission product region is worse than 
in the spallation region; therefore, development of better 
models for fission-fragment formation is of first priority. 
The new CEM2k code developed recently at Los Alamos 
p2[ motivated by the recent GSI data Q agrees best with 
our data in the spallation region, of the codes tested. But 
CEM2k is inapplicable in the fission-product region, as 
to date it has no model of fission-fragment formation. 
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